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Abstract There is accumulating evidence for mass ejection in low-mass X-ray binaries 
(LMXBs) driven by radio pulsar activity during X-ray quiescence. In this paper we con- 
sider the condition for mass ejection by comparing the radiation pressure from a millisec- 
ond pulsar, and the gas pressure at the inner Lagrange point or at the surrounding accretion 
disk. We calculate the critical spin period of the pulsar below which mass ejection is al- 
lowed. Combining with the evolution of the mass transfer rate, we present constraints on 
the orbital periods of the systems. We show that mass ejection could happen in both wide 
and compact LMXBs. It may be caused by transient accretion due to thermal instability in 
the accretion disks in the former, and irradiation-driven mass-transfer cycles in the latter. 

Key words: pulsars: general — binaries: general 
1 INTRODUCTION 

Currently there are 144 millisecond pulsars (MSPs) (pulse periods p < 10 ms) in the ATNF pulsar 
database, 86 of which are located in binary systems. MSPs are traditionally considered as the de- 
scendants of old neutron stars which were spun up by accretion from their companions in low-mass 
X-ray binaries (LMXBs) (Bhattacharya & van den Heuvel 1991 for a review). Typically, accretion of 
AM ~ Q.IMq mass is sufficient to accelerate a slowly rotating neutron star to milliseconds (Burderi et 
al. 1999). 

More recent investigations demonstrate that perhaps most of the current LMXBs have evolved 
from systems with intermediate-mass (<; 1.5M ) donor stars, i.e., intermediate-mass X-ray binaries 
(IMXBs) (Davies & Hansen 1998; King & Ritter 1999; Podsiadlowski et al. 2002; Pfahl et al. 2003). 
Because the donor star is initially more massive than the neutron star, mass transfer proceeds rapidly 
on a timescale S5 10 6 years. Little mass is accreted during this phase since the mass transfer rate is 
usually much higher than the Eddington limit for the neutron star, until the donor mass becomes com- 
parable with the neutron star mass and the binary evolves to be an LMXB (Tauris & Savonije 1999; 
Podsiadlowski et al. 2002). It is believed that in this later phase the neutron star accretes sufficient mass 
and experiences the recycling process. 

There is evidence that mass transfer in the LMXB phase is also highly nonconservative. The com- 
panion star of MSPs is usually a white dwarf of mass ~ 0.1 — 0.4M Q , the progenitor of which is 
expected to be a star with initial mass J; IMq (Tauris & Savonije 1999; Pfahl et al. 2003). So the 
lost mass from the companion is about 0.6 — O.9M . If the mass transfer is conservative, the expected 
minimum mass of the MSPs would exceed 2M@. However, measurements of the pulsar masses suggest 
that only part of the transferred mass was accreted by the neutron stars (Thorsett & Chakrabarty 1999; 
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Zhang et al. 2011: Kizilatan, Kottas, & Thorsett 2010. See also Table 1 for a list of binary MSPs with 
measured masses), except that in a few case the masses of MSPs are as high as ~ 2M (e.g., PSRs 
B1957+20, van Kerkwijk & Breton 2011; J0751+1807, Nice et al. 2005; J1614-2230, Demorest et al. 
20103). So in most of the systems, the mass transfer must have been nonconservative. 

The mechanism of mass ejection in X-ray binaries was first introduced by Illarionov & 
Sunyaev(1975) as the "propeller" effect, i.e., the accretion flow is centrifugally prohibited at the mag- 
netosphere of the neutron star, at the cost of the slow down of star's spin. As pointed out by Burderi et 
al. (2001), the efficiency of the propeller effect is at most ~ 50%, and the pulsars should still be more 
massive than observed. It was later suggested that, if the mass transfer rate varies by large amplitude, 
the neutron star may become an MSP when the accretion rate is very low, and the energetic pulsar wind 
may disrupt the accretion disk around the neutron star, so that the transferred mass will escape from the 
binary across the inner Lagrangian point L±, which we call pulsar-driven mass ejection (Burderi et al. 
2001; Burderi, D'Antona & Burgay 2002). Since this process occurs at L\, at which the binding energe 
is very small, the efficiency of ejeciton may reach unity. An alternative mechanism is that the energetic 
wind or hard X-ray radiation from the MSP may evaporate it's low-mass companion (van den Heuvel 
& van Paradijs 1988; Ruderman, Shaham & Tavani 1989; Shaham & Tavani 1990; Shaham & Tavani 
1991;Podsiadlowski 1991). 

Direct observational evidence for mass ejection induced by pulsar activity may come from the 
"black widow pulsars". They are MSPs in binary systems, undergoing very wide eclipse, implying ob- 
scuration by intense wind from the secondary (Fruchter, Stinebring, & Taylor 1988). Another interesting 
example is PSR J1740— 5340, an eclipsing MSP with a spin period of 3.65 ms and orbital period of 32.5 
hr, located in the globular cluster NGC 6397 (D' Amico et al. 2001). Long lasting and sometimes irreg- 
ular radio eclipses, and the shape of the optical light curve demonstrated the presence of matter around 
the system (Ferraro et al. 2001), suggesting that PSR J1740— 5340 is an example of a system in the 
pulsar-induced ejection phase. 

In this paper we investigated the conditions for mass ejection in LMXBs caused by the turn-on of 
pulsar activity. The basic idea is similar to that in Burderi et al. (2001). However, we consider direct mass 
ejection at both L\ and L2, which requires that the pulsar's radiation pressure is strong enough to drive 
material out of these points, and was ignored in Burderi et al. (2001). For mass ejection due to disruption 
of the accretion disks by the pulsar's pressure, we adopt a more realistic accretion disk model to evaluate 
the disk pressure. In section 2 we calculated the critical periods p CI at which the pulsar wind pressure 
P psr equals the donor's gas pressure Pli at the L\ point, and Pdisk in the accretion disk respectively. 
The condition for outflow from L2 is also considered. In section 3 we put possible constraints on the 
orbital periods by combining the evolution of LMXBs. We briefly discuss the implications of our results 
in section 4. 

2 THE CRITICAL PERIODS 

2.1 Mass Flow Through the Inner Lagrange Point L\ 

We consider a semi-detached binary containing a neutron star and a donor star that fills it's Roche lobe. 
The atmospheric material is overflowing through L\ towards the neutron star. At L\ the gas flow is 
confined in a "nozzle" with radius H and moves with a speed close to the sound velocity (Frank, King, 
& Raine 2002). The mass transfer rate is expressed as 

M ~ TrH 2 c s (Lx)p(Li), (1) 

where c s (Li) and p(Li) are the local sound velocity and the gas density at the L\ point. To estimate the 
magnitude of H, we consider the balance between the Roche potential and the mean kinetic energy of 



In a recent paper, Tauris et al. (201 1 ) suggested that PSR J1614 — 2230 must have been born with a mass significantly exceeds 
1.4 M Q . 
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the material, and find that matter escapes in a patch of radius (see also Frank, King, & Raine 2002; Li 
et al. 2010) 

c.(Li) 



H 



where 



(1 + ? )(0.5 - 0.2271og 10 qf (1 + g)(0.5 + 0.227 log 10 qf 



(2) 



(3) 



where q is the mass ratio of the donor and the neutron star, ranging from ~ 0.1 to ~ 3 for I/LMXBs. 
The mass transfer rate can then be written as 

M = 1.03 x 10 e P^ h hl 4(L 1 )p(L 1 )f^ 1 (q)g S ~\ (4) 

where P or b,hr is the orbit period of the system in units of hour. Consequently, the gas pressure at L\ is 

P L1 = = 6.15 x 10 3 M_ 10 P-" hr c s - 6 1 /i dyncm" 2 , (5) 

where M- W = M/10 -10 Af yr _1 , andc s6 = Cg/10 6 cms -1 . 



2.2 Equilibrium at L\ and L 2 

Once the material from the donor star is captured by the neutron star at the circularization radius, an 
accretion disk will form due to dissipative processes and internal torque in it. Sufficiently long time 
accretion on to the neutron star can spin it up to millisecond periods. If the accretion rate drops due 
to some reasons, the inner radius of the accretion disk may move outside of the light cylinder (with 
radius Rlc) of the neutron star, and the neutron star changes from an accreting X-ray source to be 
a rotation-powered pulsar, emitting out radio, optical, X-ray, 7-ray photons and high energy particles. 
If this process is dominated by the magnetic dipole radiation, the pulsar wind pressure at L\ can be 
expressed as 

P psr (i!) = 8.23 x 10 2 B 2 s8 p-^lp^R e n6 mf /3 f^ 2 (q) dyncm" 2 , (6) 

where mi and R n § are the mass and radius of the pulsar in units of solar mass and 10 6 cm respectively, 
B s s the surface magnetic field in units of 10 8 G, p ms the pulsar period in units of millisecond, and 

f 2 (q) = (1 + g)V3(0.5 - 0.2271og 10 q). (7) 

When the pulsar wind pressure is larger than the gas pressure (Pli) at L\, the accretion process is 
prohibited, and mass ejection will happen. Combining Eqs. (5) and (6) we can get a critical period at 
which material flow through L\ will be ejected, 

r> ~ fn fil ™ e \ f- 1 /4 f -l/2,V-l/4 l/4 R l/2 p l/6 R 3/2 -1/6 m 

p cr ~ (O.ol m&j /j J 2 M_ w c s6 B sS P orb;hr -K n6 m . (8) 

When the secondary star evolves to be of very low mass, it may be the outer Lagrangian point L2 
rather L\ with the lowest potential, due to the pulsar's radiation pressure force on the Roche potential of 
the secondary. In this case, Roche-lobe overflow will not occur towards the pulsar, but into a circumbi- 
nary disk/outflow (as is indeed observed for PSR 1957+20). According to Phillips & Podsiadlowski 
(2002), the condition for outflow from L 2 can be defined as 8 max = qL/LEdd becomes larger than a 
critical value S cr a (where L and L Edd are the irradiating pulsar's luminosity and the Eddington lumi- 
nosity of the secondary, respectively), or equivalently 
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where 



L 



2Rl 



B'i I — 



3c 3 J " \ p 

and k is the mean photospheric opacity. The resulting critical spin period is 



(10) 



Pcr,L2 = 27T 



2K V /4 



/ kB 2 

\4ircGM5 crit 3c 3 J 



(ID 



Pcr,L2 - (0.48 ms) K 



1/4 
0.4 7 



-l/4 R l/2 R 3/2 



(12) 



with S cr it = 0.0564 when q = 0.1 (Phillips & Podsiadlowski 2002). Here K0.4 = k/0.4 cm 2 g _1 may 
be much larger than unity for the photosphere of very low-mass stars, so Eq. (12) presents a lower limit 
of the critical period. Obviously, depending on the mass ratio, this criterion may be more relevant for 
very low-mass binaries. 

2.3 Equilibrium with the disc 

Now we consider possible disruption of the accretion disk by the pulsar's pressure. The standard geo- 
metrically thin a-disk model was established by Shakura & Sunyaev (1973). However, there is obser- 
vational evidence showing deviation from the theoretical predictions - in cataclysmic variables (CVs) 
and LMXBs the outer disk edge seems to be thick and structured (Shafter & Misselt 2006; Hakala et 
al. 1999). Begelman & Pringle (2007) argued that, the magnetorotational instability (MRI), which suc- 
cessfully accounts for angular momentum transport in accretion disks (Balbus & Hawley 1998), can 
amplify the toroidal magnetic field to a point at which magnetic pressure far exceeds the combined gas 
and radiation pressure in the disk. This additional pressure support makes the disk thicker than in the a- 
disk model, in accordance with observations. In this magnetic-dominated disk model, the Alfven speed 
associated with the toroidal field is roughly the geometric mean between the Keplerian speed and the 
sound speed. Thus the magnetic pressure is Pg ~ pc g Vk- The total disc pressure can be derived to be 
(see appendix) 



P, 



disk 



P B = 1.02 x 10 13 a 



17/18 m 61/36 fl>8/9 p-91/36 R -l 



R 



n 6 



Rn6 



1/2 



-1/9 



dyn cm 



(13) 

where R$ is the disk radius in unit of 10 6 cm. As pointed out by Burderi et al. (2001), if the spun-up 
pulsar is switched on, the accretion disk may be truncated at a radius i? s top due to the energetic pulsar 
wind. Once i? s top is smaller than the outer radius i? ou t of the accretion disk, the overflowed material 
through L\ will be ejected, i.e., the system will enter the mass ejection phase. 

Traditionally the outer radius of the disk is approximated to be ~ 70% — 90% of the Roche lobe 
radius Rn of the primary (Frank, King, & Raine 2002). Here we take i? ut — O.SRli- Combining 
Eq. (6) and (13) , the critical pulsar spin period for mass ejection is 



{IM^a^m^R^B^MZl^P^P. 



(14) 



where 



h = (1 + 3) 



19/432 



1 - 0.462 



1/3 



91/144 



(0.5- 0.227 log 10 q) 



-1/2 
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2.4 Accretion Equilibrium Period 

In the above subsections we show that mass ejection can take place only when the spin period of the 
neutron star is shorter than either p CI or p' CT . The neutron star's spin is accelerated by both mass accre- 
tion and magnetic field-disk interaction during the mass transfer process, until an equilibrium period 
is reached, at which there is no net torque exerted on the neutron star. According to Ghosh & Lamb 
(1979), the spin-up/down rate depends on the fastness parameter u> s = CI s /CIk(i"o)> i- e -> the ratio of the 
spin angular velocity of the neutron star and the Keplerian velocity of plasma at the inner disk radius 
ro. When ui s approaches the critical fastness parameter lj c , the torque is zero and the spin-up process 
ceases. Assume ro = 0.5i?A, where Ra is the traditional Alfven radius for spherical accretion 

B 2 s Rl \ 2/7 
AIV2GAIJ 

the equilibrium spin period is 

Pcq ~ (1.95 ms) u^R^Blgm-^MZll 7 , (15) 

where ujqj = ui c /0.7. The condition p cq < p CI or p' cr will constrain the systemic parameters that allow 
mass ejection to happen in LMXBs. 

The neutron star will reach its equilibrium only when the mass transfer rate keeps constant for 
sufficient long time. When the mass transfer rate changes, the spin period will evolve to new equilibrium 
on a time scale t sp in ~ 2ttIv s /AI(GMR) 1 / 2 ~ 2.3x lQpp^MZio, wnere H s the momentum of inertia 
of the pulsar. Roughly speaking, if the timescale of the change of the mass transfer rate > i sp i n , the 
spin period of the neutron star is always close to the equilibrium period with current M. If M varies 
with tj^ <?C tspin, the spin period remains to be close to the previous equilibrium value. 




3 BINARY MSP SYSTEMS THAT FAVOR MATERIAL EJECTION 
3.1 Mass Transfer Mechanisms 

In this section we investigate the possible influence of the mass transfer rate. As mentioned above, 
recent studies showed that most LMXBs are likely to originate from IMXBs. During the initial TMXB 
to LMXB' phase, the mass transfer process proceeds on a (sub)thermal timescale, the average mass 
transfer rate is very high (often super-Eddington for relatively massive companions), so the mass transfer 
is highly unconservative and the accreted mass by the neutron star is very small. Thus the recycling 
process should mainly occur in the later LMXB phase. From Eqs. (8), (12) and (14) one can see that 
mass ejection due to pulsar activity can take place only when the neutron star's spin period is less 
than a few milliseconds. This is difficult to achieve except for LMXBs and close IMXBs (with case A 
mass transfer); the outcome of wide IMXB evolution is mildly recycled MSPs with CO WD companions 
(Tauris 201 1). In the former case the LMXB evolution from IMXBs is similar to that of original LMXBs 
(Podsiadlowski et al. 2002). For these reasons we only consider the evolution of mass transfer rate in 
LMXBs. 

It's well known that mass transfer in LMXBs is driven by nuclear expansion of the donor and 
orbital angular momentum loss due to magnetic braking and gravitational radiation. By comparing the 
nuclear evolution time of the donor t nu and the angular momentum loss timescale t am \, we can divide 
the evolutionary tracks of LMXBs into the following cases (King, Kolb, & Burderi 1996). 

Case 1 t nu <C t a mi: In this case the donor has evolved off the main sequence before angular momen- 
tum loss shrinks the orbit sufficiently to cause mass transfer. So the donor must be a giant or subgiant, 
with mass of > 0.8M Q in order to evolve off main sequence within Hubble time. It's known that the 
structure of a low-mass giant is determined by the mass of its helium core (Webbink et al. 1983). The 
mass transfer rate can be expressed by the following relation (King 1988), 

m_ 10 c o.2ii£»3x- 4 y- 47 . (16) 
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Fig. 1 The relation between P rb,hr and the mass ratio q for LMXBs in case 1 evolution when 
mass ejection occurs. The solid and dashed lines are for p CT = p cq , and p' CI = p cq . Above 
the curves the mass ejection may happen. The dotted line represents the condition when the 
donor fills its Roche lobe on the main sequence. Here and in the following we adopt m = 1.4, 
R n6 = 1, B b8 = 5 and a = 0.1. 
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Fig. 2 The relation between critical spin period and the mass ratio q for LMXBs in case 1 
evolution. The black, blue and red lines describe p CT , p' CI , and p cq , respectively. The dashed and 
solid lines are for orbital period P rb,hr = 100, 1000, respectively. 

Substitute Eq. (16) into Eqs. (8), (11), and (15), when p cr = p cq and p' CI = p cq we obtain the minimum 
of the orbital period for mass ejection 

P — 7K Qfi , ,- 2 - 94 ™- 2 - 4 0„-0.78 R 1.08 f 0.76 f 1.50.-0.76 p3.21 n 7 x 

r'oibM — 75.9b u 7 m q B sS } 1 } 2 c s6 H n6 , (17) 

and 

p/ _ o 00 -3.45„-0.83„-1.07 p2.83 R 1.24_.-3.78 f -3A5 n o\ 

Figure 1 shows the minimum P rb,hr and P^ rh hr against the mass ratio q in the solid and dashed lines, 
respectively. Note that the minimum orbit period that allows direct mass ejection at L\ is ~ 5 days, 
while mass ejection with disk disruption can occur in relatively compact systems. The final product of 
this kind of evolution is a wide binary with an MSP and a low-mass helium white dwarf. 

We also plot p cr , p' CI and p cq for constant orbital period with black, blue and red lines in Fig. 2 as a 
function of mass ratio q. Note that in the cases of P rb,hr = 100 (dashed lines) and 1000 (solid lines), 
p' CI is always larger than p cr . 

Case 2 t nu ~ t aml and i nu ^> t am i: In these situations the mass tranfer is driven by angular momen- 
tum loss of the system, and 

Ma J 
M 2 ~ J' 
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L°g 10 p o* (hour) 

Fig. 3 This figure compares p CI , p' c[ , and p eq as a function of orbital period in case 2 mass 
transfer. Here p CT , p' cr , and p eq are plotted in black, blue, and red lines, respectively. The dashed 
and dotted lines are for q — 0.5 and 0.1, respectively. The dot-dashed line represents p CI ^2 for 
q = 0.1. 

where J and J are the orbital angular momentum and its derivative, respectively. In the case of t nu ~ 
iami the donor is mildly evolved when it fills it's Roche lobe, and the radius is slightly bigger than during 
the main sequence. The orbit period of the system will not change much during the whole evolution. 
When t nu t am \, the angular momentum loss of the system is so rapid that mass transfer occurs when 
the secondary is still on the main sequence, and further nuclear evolution of the secondary is frozen. 
In this condition the binary orbit will shrink due to the angular momentum loss by magnetic breaking 
and gravitational radiation. Hence the orbital periods are generally smaller than in case 1 where the 
secondary is a (sub)giant. 

Use the standard forms of gravitational radiation losses (Landau & Lifschitz 1958) and of magnetic 
braking (Verbunt & Zwaan 1981), the mass transfer rate can be written as (King et al. 1996), 

M_ 10 = 552.81P OT ^ h y/ 3 m 5 / 3 + 130.65P OT ^ h V m 8 / 3 , (19) 

where the first and second terms on the right hand side are for magnetic braking and gravitational 
radiation, respectively. Figure 3 shows p CI , p' CI , and p eq as a function of P or b in black, blue, and red 
lines, with the mass transfer rate given by Eq. (19). The dashed and dotted lines are for q = 0.5 and 0.1 
respectively. In the case of q = 0.1 we use Eq. (19) with donor mass = 0.3 Mq to calculate p oq , but 
remove the first term (i.e., corresponding to magnetic braking) in calculating p CT and p' cr . The reason is 
that magnetic braking is assumed to vanish abruptly when the secondary mass decreases to ~ 0.3 M© 
(i.e., becomes fully convective), and the subsequent evolution is driven by gravitational radiation only. 
However, it will take a long time (<; 10 9 yr) for the pulsar to reach the new equilibrium period because 
of the reduced mass transfer rate. So the spin period remains close to the previous one. We have set the 
limitation on the orbital period 0.11P or b,hr > m 2 to guarantee that the secondary fills its Roche lobe 
before it evolves to be a subgiant. Figure 3 shows that p CT & p cq < p' CI , which means that mass ejection 
via disk disruption might occur. But a firm conclusion is difficult to reach, because of the uncertainties in 
estimating the critical periods. By the way, outflow from L2 seems to require spin periods considerably 
lower than p cq for q = 0.1, disfavoring this possibility for LMXBs. However, this scenario might work 
in binaries with extreme mass ratio (q <C 1, like PSR B1957+20), in which the value of <5 cr ,t is much 
smaller (Phillips & Podsiadlowski 2002). 

Note that p cq < p cx or p' cx is the necessary condition for mass ejection, while switch-on of pulsar 
emission is critical in expelling the matter overflowing the Roche lobe. This occurs with a temporary re- 
duction of the mass accretion rate, so that the accreting plasma moves out beyond the light cylinder of the 
neutron star, and the neutron star becomes a generator of magneto-dipole radiation and relativistic par- 
ticles. The thermal disk instability model (Lasota 2001) is usually invoked to trigger the X-ray outbursts 
in LMXBs, which may cause the switch-on of the radio pulsar during X-ray quiescence. According to 
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Table 1 Properties of binary MSPs with measured masses 



PSR Name 


P (ms) 


-Porb.hr (hr) 


M P (M ) 


M 2 (M Q ) 


Ref. 


B1516+02B 


7.95 


164.6 


2.08 ±0.19 


0.13* 


1 


B 1855+09 


5.36 


295.2 


1.6 ± 0.2 


o.258tg:E 


2 


B 1957+20 


1.61 


9.17 


2.40 ± 0.12 


0.035±S;£ 


3 


J0024-7204H 


3.21 


60.9 


1 41+0.04 

— 08 


0-18±g-gg 


4 


J0437-4715 


5.76 


137.8 


1.76 ±0.20 


0.254 ±6.018 


5 


J075 1+1807 


3.48 


6.3 


1.26 ±0.14 


0.19 ±0.03 


6 


J1012+5307 


5.26 


14.5 


1.68 ±0.22 


0.16 ±0.02 


7 


J1614-2230 


3.15 


208.8 


1.97 ±0.04 


0.4* 


8 


J 17 13+0747 


4.57 


1627.8 


-. ro+0.08 
- L -' J " J — 0.06 


0.33 ± 0.04 


9 


J1738+0333 


5.85 


8.52 


1 55 ± 55 


2 ± 05 


10 


T 1740- 5 340 A 


3.65 


32.5 


1 53 ± 1 9 


0.22* 


1 1 


J 1903+0327 


2.15 


2284.2 


1.74 ± 0.04 


1 05 ± 02 


12 


J 1909-3744 


2.95 


36.7 


-, 47 +o.03 

J-- 1 *' -0,02 


0.2038 ± 0.0022 


13 


J1911-5958A 


3.27 


20.6 


1 4+O.I6 


0.18 ±0.02 


14 


J1748-2446I 


9.57 


31.87 


1 87+ ' 32 

1 -°'-0.07 


0.24* 


15 


J 1045-4509 


7.47 


97.92 


< 1.48 


0.13 


16 


J 1804-27 18 


9.34 


266.4 


< 1.73 


0.2 


16 


J20 19+2425 


3.93 


1836 


< 1.51 


0.32 - 0.35 


17 


J 1023+003 8 


1.69 


4.75 


1.0-3.0 


0.14 - 0.42 


18 


J1738+0333 


5.85 


8.50 


1.6 ±0.2 


0.2 


19 


J0024-7204I 


3.49 


5.52 


1.44 


0.15 


20 


J1518+0204B 


7.95 


164.64 


2.08 ±0.19 


> 0.13 


21 


J1824-2452C 


4.16 


193.87 


< 1.367 


> 0.26 


21 


J0514-4002A 


4.99 


450.96 


< 1.52 


> 0.96 


22 


J 1748-202 11 


9.57 


31.87 


1.3 ±0.02 


0.24 


15 



Notes: *: Medium companion mass from ATNF pulsar database. l.Freire et al 2008a, 2.Splaver 2004, 3. van 
Kerkwijk et al. 2011, 4.Freire et al. 2003, Manchester et al. 1991, 5.Verbiest et al. 2008, 6.Nice et al. 2008, 
7.Lange et al. 2001, 8.Demorest et al. 2010, 9.Splaver et al. 2005, lO.Freire et al. 2008b, 1 l.Kaluzny et al. 2003, 
12.Champion et al. 2008, 13 Jacoby et al. 2005, Hotan et al. 2006, 14.Bassa et al. 2006, 15.Ransom et al. 2005, 
16.Thorsett & Chakrabarty 1999, 17.Nice, Splaver, & Stairs 2001, 18.Archibald et al. 2009, 19.Jacoby 2004, 
20.Manchester et al. 1991, 21. Zhang et al. 2011, 22.Freire et al. 2007 

King et al. (1996), long-period LMXBs are likely to be transient; for LMXBs with P orb ^ 2 days, 
most neutron star systems will be persistent X-ray sources unless they have extreme mass ratios (i. e., 
q < 0. 18Porb.hr)- Consider the calculated results in Figs. 1-3, we conclude that mass ejection is more 
likely to occur in wide systems with light secondary (or in the late the evolution of long-period LMXBs) 
via disruption of the accretion disk inside the Roche lobe of the primary. 

4 SUMMARY AND DISCUSSION 

We can summarize the above results for the conditions of pulsar-driven mass ejection as follows. 

(1) Pulsar-driven mass ejection is likely to occur in wide LMXBs (in case 1 evolution), in which the 
thermal instability in the accretion disk can result in large variations in the accretion rate. 

(2) LMXBs in case 2 evolution may barely be subject to mass ejection. 

(3) If there is mass ejection, it may be caused by the disruption of the accretion disk leading to 
outflow through the L\ point. When the secondary is of extremely low mass, outflow from the point 
is also possible. 

An important factor that was neglected in the above sections is irradiation in LMXBs caused by 
accretion-generated X-rays. It was realized that irradiation of the donor star or the accretion disk can not 
only change the optical appearance of LMXBs (King & Ritter 1999), but also their outburst properties 
(van Paradijs 1996), and possibly the long-term evolution of the donor (Podsiadlowski 1991). Especially 
irradiation of the donor star can destabilize the mass transfer, and lead to mass transfer cycles (Hameury 
et al. 1993), which differs drastically from the evolution we considered above: mass transfer is spas- 
modic with phases of high mass transfer driven by the thermal expansion of the convective envelope 
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of the irradiated donor alternating with phases with low or no mass transfer, during which the donor 
readjusts towards thermal equilibrium of the un-irradiated star (Biining & Ritter 2004). As pointed out 
by Ritter (2008), the effect of irradiation may be important in compact rather long-period LMXBs, since 
in the later the irradiation resulting from accretion is intermittent due to disk instability. However, even 
in the former systems, the details of how irradiation of the donor influences the mass transfer process are 
very complicated, and have not been included in the calculations of the secular evolutions of I/LMXBs 
in a self-consistent way (Pfahl et al. 2003). It is then difficult to estimate the efficiency of mass ejection 
due to irradiation. On one hand, during the short "high" state the mass transfer rate is enhanced and 
larger than the secular one, leading to shorter p eq , which is in favor of mass ejection (In Fig. 3, if the 
mass transfer rate is increased by a factor of 100, both p cq < p cr and p cq < p' cr will be satisfied). On 
the other hand, during the long "low" state, the mass transfer rate becomes much lower than the secular 
one, so that less mass can be blown off. This is in contrast with the limit cycles due to disk instability, 
in which the accretion rate varies by a large factor but the mass transfer rate remains nearly unchanged. 

The current masses of MSPs may reveal possible evidence of mass accretion and ejection during the 
previous LMXB evolution. In Table 1 we list the parameters of binary MSPs with measured masseo In 
general, long-period (P rb 20 — 30 days) MSPs have masses around 1.4 M , suggesting that large 
amount of the transferred mass was lost. In narrow systems (P rb < 2 — 3 days) the pulsar masses 
distribute from ~ 1.3 Mq to ~ 2.4 Mq, indicating that both efficient mass accretion and mass ejection 
are possible, depending on the properties of the individual sources. However, a large fraction of them 
also have masses not far from 1.4 Mq. We speculate that irradiation-driven mass transfer cycles may 
help drive off the transferred mass from their companions in these systems. 

Besides the standard recycling scenario, another possible way for mass ejection in short-period 
LMXBs is companion exchange. This should happen only in the environment with high stellar density 
like globular clusters. The incidence of black widow pulsars in globular clusters is known to be far 
higher than in the field. This leads King et al. (2003) to suggest that the MSPs in dead wide binaries 
with white dwarf companions have exchanged them for normal stars. Encounters and tides bring these 
new companions into tight orbits. Due to intense accretion during the first mass transfer phase, the 
neutron star's spin period may satisfy the condition p cq < p cr or p' CI with current mass transfer rate, 
resulting in mass ejection from the binary. 

Finally, the condition for mass ejection may be satisfied if the accreting star in LMXBs is a strange 
star rather a neutron star. In the former case the spin periods could be submilliseconds (e.g., Frieman & 
Olinto 1989; Gourgoulhon et al. 1999; Possenti et al. 1999), and less than the critical periods given by 
Eqs. (8) and (14). These exotic objects might be formed by accretion-induced collapse of white dwarfs 
in binary systems (Du et al. 2009). 
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Appendix A: CALCULATION OF THE DISK PRESSURE 

In the model of magnetic dominant disk Begelman and Pringle argued that the amplified toroidal mag- 
netic field pressure is ~ pc g Vk- They calculated the density inside the disc is 

£ 2 c 2 / c\ 3/2 M _! / i? V 9/4 

-a ( — ) (A.l) 

here k is the opacity, Msdd = AivGM/kc is the Eddington accretion rate and R g = GMjc?. In this 
model the viscosity and scale height of the disc is given by v = aHvA and H/R ~ (cg/vk) 1 ^ 2 instead 



2H 3GMk\cJ M E dd \R 



2 PSR J1903+0327 is a peculiar MSP with a main-sequence companion star, and is not considered here. 
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of v = aHcg and H/ R = {c? g + cf.) 1 I 2 /vk in SS model. After these replacement then use the equation 
of angular momentum 



v M 



and energy 



4a r 
3r 



-T,' 



we can calculate the disk pressure is 

17/18 (GM)6V3« f 3k 



6lTC 2 V 327T 2 (T 



3GMM 
8ttR 3 



-Vis M 



1/2' 



M 8 / 9 / k 



R n \ jirrip 



1/2 



(A.2) 



(A.3) 



-2/9 



^-91/36 



J2. 



1/2' 



-, -1/9 



(A.4) 
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